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ABSTRACT 

Detailed  interferometric  measurements  oi  the  flow 
near  the  leading  edge  of  an  oscillating  airfoil  offer 
the  first  detailed  experimental  quantification  of  the 
locally  compressible  flow  field  that  surrounds  an  os¬ 
cillating  airfoil  at  moderate  subsonic  Mach  numbers. 
Interferograms  obtained  zby  a  specially  adapted  real¬ 
time  point-diffraction  interferometry  technique  have 
revealed  significant  characteristics  of  this  complex, 
and  very  rapidly  varying,  locally  supersonic  flow.  In¬ 
stantaneous  pressure  distributions  determined  from 
these  interferograms  document  the  effect  of  unsteadi¬ 
ness  on  the  leading-edge  flow  environment. 

Nomenclature 

Cp  pressure  coefficient 

c  airfoil  chord 

f  frequency  of  oscillation,  Hz 

k  reduced  frequency  = 

L  test  section  span 

Moo  free  stream  Mach  number 

n  refractive  index 

nr  refractive  index  at  reference  conditions 

no  refractive  index  at  atmospheric  conditions 

n  distance  along  airfoil  surface 

Uoo  free  stream  velocity 

X  chordwise  distance 

Q  angle  of  attack 

APL  average  path  length  difference 

f  fringe  number 
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Ao  wavelength  of  the  laser 

p  density 

Po  density  at  atmospheric  conditions 

pr  density  at  reference  conditions 

<i>  phase  angle  of  oscillation 

w  circular  frequency,  radians/sec 

1.  INTRODUCTION 

Quantification  of  the  complex  unsteady  flow  field 
associated  with  dynamic  stall  continues  to  be  of  sig¬ 
nificant  importance.  The  present  inability  to  con¬ 
trol  and/or  modify  the  dynamic  stall  phenomenon 
has  limited  the  speed  and  agility  of  helicopters,  and 
has  seriously  restricted  the  use  of  unsteady  lift  for 
improvement  of  aircraft  agility  and  maneuverability. 
Many  studies  of  this  aerodynamic  phenomenon  have 
been  performed;  see  Carr^  for  summary.  Recent  stud¬ 
ies  using  the  facility  described  in  the  present  paper 
have  shown  that  compressiblity  effects  can  and  do 
play  a  significant  role  in  the  development  of  the  un¬ 
steady  flow  field  on  rapidly  pitching  airfoils^’^.  How¬ 
ever,  truly  quantitative,  simultaneous  documentation 
of  the  compressible  flow  field  on  and  away  from  the 
surface  of  the  airfoil  has  not  been  available  up  to  the 
present.  Surface  measurements  have  been  made  in  the 
past'*’^,  but  knowledge  of  the  flow  behavior  away  from 
the  surface  is  also  necessary  before  control  and/or 
modification  of  the  stall  process  can  be  implemented. 
The  work  presented  in  this  paper  supplies  some  of  this 
critical  information. 

2.  FACILITY  AND  INSTRUMENTATION 

The  present  study  was  performed  in  the  Com¬ 
pressible  Dynamic  Stall  Facility(CDSF)  at  the  NASA 
Ames  Research  Center  Fluid  Mechanics  Laboratory. 
The  CDSF  is  specifically  designed  for  study  of  dy¬ 
namic  stall  over  a  range  of  Mach  numbers,  using  non- 
intrusive  optical  flow  diagnostic  techniques.  It  is  op¬ 
erated  as  a  part  of  the  in-draft  tunnel  complex  at  the 
FML(for  details  see  Carr  and  Chandrasekhara,  Ref¬ 
erence  6).  The  CDSF  is  unique  in  that  the  airfoil  is 
supported  between  two  2.54  cm  thick  optical  quality 
glass  windows  by  pins  that  are  smaller  than  the  local 
airfoil  thickness.  Thus,  the  entire  flow  field  includ- 
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ing  the  airfoil  surface  can  be  viewed  unobstructed  by 
any  support  nieclianism(Figure  1).  This  permits  the 
study  of  the  flow  at  the  surface  near  the  leading  edge, 
where  the  dynamic  stall  vortex  forms,  as  well  as  the 
flow  field  away  from  the  airfoil.  Details  of  the  research 
findings  from  the  earlier  work  in  this  facility  can  be 
found  in  References  2,  3  and  6.  The  present  tests  were 
performed  for  a  =  10"  -  10°sinu;t,  0.28  <  Mco  <  0.47, 
and  0.0  <  k  <  0.10. 

3.  FLOW  VISUALIZATION 

3.1  Point  Diffraction  Interferom  try  Technique 

Poiiit  diffrci^tiuii  interfcrometry(PDI)  is  a  tech¬ 
nique  which  permits  recording  of  real-time  interfero- 
grams  of  a  flow'"*'’.  The  PDI  technique  utilizes  the 
ability  of  a  point  discontinuity(in  the  form  of  a  pinhole 
or  an  opaque  spot),  located  at  the  image  of  a  point 
sourcc( Fourier  transform  plane),  to  diffract  a  portion 
of  the  incident  light  into  a  spherical  reference  wave 
front.  Light  passing  through  the  ph<ise  object  under 
test  is  focused  onto  a  sheet  of  partially  transmitting 
material  containing  the  pinhole(see  Figure  2).  A  por¬ 
tion  of  this  light  is  diffracted  by  the  pinhole  creating 
a  reference  wavefront.  The  resulting  reference  wave- 
front  then  interferes  with  the  undiffracted  light  that 
passes  around  the  pinhole.  The  interference  fringes 
which  are  produced  by  this  method  provide  a  contour 
of  optical  path  difference  caused  by  density  variations 
in  the  flow  field.  Since  the  mathematical  derivations 
associated  with  PDI  are  well  documented  in  the  refer¬ 
enced  papers'"®,  the  present  paper  will  only  address 
the  physical  implementation. 

3.2  Implementation 

In  the  present  application,  the  primary  optics  of 
an  existing  schlieren  system  are  used,  with  a  pulsed 
Nd-YAG  laser  replacing  the  conventional  spark  as  the 
light  source,  and  the  above  mentioned  point  diffrac- 
tor  replacing  the  usual  knife  edge.  Figure  3  shows 
a  schematic  of  the  arrangement.  Laser  light  is  ex¬ 
panded  to  fill  a  portion  of  the  first  schlieren  mirror, 
where  it  reflects  through  the  test  section  and  finally 
is  refocused  onto  the  diffracting  plate  by  the  second 
schlieren  mirror.  A  camera  placed  at  the  image  plane 
of  the  flow  field  records  the  interference  pattern. 

Diffracting  spots  and  pinholes  were  made  by  plac¬ 
ing  a  photographic  emulsion  precisely  at  the  focus 
of  the  second  schlieren  mirror.  With  no  flow  in  the 
test  section,  laser  energy  was  incrccised  until  a  spot 
was  charred  on  the  emulsion,  or  until  a  pinhole  was 
burned  through  the  emulsion.  Although  some  of  the 
earlier  data  was  obtained  with  diffracting  spots,  pin¬ 
holes  were  used  in  the  present  experiment  since  the 
interferograms  produced  using  pinholes  were  brighter 
and  of  superior  quality.  Once  the  pinhole  was  made, 
it  was  left  in  place,  and  the  laser  energy  decreased 
to  a  level  thf-t  adequately  exposed  the  camera  film 
without  destroying  the  pinhole. 

The  photographic  plate(Agfa  8E75)  used  to  make 
the  diffracting  pinhole  was  processed  with  a  weak  de¬ 


veloper  and  then  exposed  uniformly  with  light  in  or¬ 
der  to  make  the  emulsion  darker(plates  with  emul¬ 
sions  darkened  to  optical  densities  of  approximately 
one  were  used  for  the  flow  situations  encountered  in 
this  test);  no  other  processing  was  required.  The  in¬ 
creased  absorption  is  necessary  to  attenuate  the  light 
passing  around  the  pinholefobject  wave)  to  a  level 
equal  in  intensity  to  the  light  which  diffracts  from  the 
pinhole  (reference  wave!,  thereby  maximizing  fringe 
contrast.  The  optical  density  necessary  depends  on 
how  much  the  light  is  disturbed  as  it  passes  through 
the  test  section.  An  increase  in  disturbance(caused  by 
stronger  density  gradients  such  as  those  that  appear 
at  high  angle  of  attack)  results  in  an  increase  in  the 
diameter  of  the  focused  beam  at  the  pinhole,  thcicLy 
decreasing  the  light  diffracting  through  the  pinhole, 
thus  requiring  a  higher  optical  density  for  the  PDI 
plate  for  optimum  fringe  contrast. 

An  alternate  method  of  matching  the  beam  inten¬ 
sities  was  used  for  certain  conditions.  This  method 
used  carefully  controlled  laser  energy  to  burn  vary¬ 
ing  sizes  of  pinholes.  Since  a  wide  range  of  density 
gradients  were  encountered  in  the  flow  fields  of  inter¬ 
est,  frequent  variation  of  either  emulsion  density  or 
pinhole  weis  required.  This  latter  method  of  match¬ 
ing  beam  intensities  required  the  least  effort,  and  was 
used  for  most  of  the  interferograms  presented. 

4.  FLOW  CHARACTERISTICS 
4.1  General  Flow  Characteristics 

The  interferogram  presented  in  Figure  4  shows  in 
detail  the  developing  compressible  flow  field  around 
the  leading  edge  of  the  pitching  airfoil  just  prior  to 
the  initiation  of  dynamic  stall.  This  figure  presents 
the  flow  at  a  =  10.0"  for  Moo  =  0.33,  and  k  =  0.075, 
and  shows  several  characteristics  of  the  unsteady  flow 
field  that  are  of  significant  interest.  First,  the  stagna¬ 
tion  point  on  the  lower  surface  is  eeisily  discernable,  at 
the  center  of  the  smallest  closed  contour  of  constant 
density,  since  the  flow  speed  increases  from  the  stag¬ 
nation  point  in  both  directions  on  the  lower  surface  of 
the  airfoil,  it  is  clear  that  the  stagnation  point  is  cen¬ 
tered  in  this  fringe.  It  should  also  be  noted  that  this 
stagnation  point  quantifies  the  density  field  for  the 
complete  interferogram,  since  the  stagnation  pressure 
of  the  surrounding  flow  is  attained  at  the  stagnation 
point  on  the  airfoil.  This  allows  quantification  of  the 
fringe-related  Mach  contours  for  the  complete  flow 
field  around  the  airfoil. 

The  flow  accelerates  eis  it  progresses  around  the 
leading  edge,  eis  can  be  seen  in  the  closely  spaced 
fringes  emanating  from  the  leading  edge  region.  A 
small  region  where  no  image  is  present  can  be  ob¬ 
served  here  at  the  leading  edge;  this  region  is  dark  due 
to  the  fact  that  much  of  the  light  has  been  so  severely 
deflected  by  the  flow  gradients  that  it  is  not  able  to 
reach  the  photographic  plate,  thus  leaving  a  dark  re¬ 
gion  on  the  film.  Considerable  effort  was  expended 
during  the  alignment  phEise  to  minimize  this  dark  re¬ 
gion;  in  addition,  during  analysis,  the  imaf^e  in  this 
region  was  further  improved  using  a  digital  enhance¬ 
ment  technique  described  below.  As  the  flow  contin¬ 
ues  to  move  around  the  airfoil,  it  continues  to  accel- 
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erate,  until  the  surface  curvature  is  reduced,  where  a 
plateau  appears  in  the  density  contours.  This  region 
is  of  significant  interest,  since  surface  pressure  mea¬ 
surements  performed  in  earlier  experiments^  suggest 
locally  supersonic  flow  in  this  region  based  on  isen- 
tropic  expansion  to  the  surface  pressure  measured  by 
that  experiment.  However,  the  constant-density  con¬ 
tours  obtained  in  the  present  experiment  show  that 
for  the  present  experimental  conditions,  this  super¬ 
sonic  region,  if  present  at  all,  is  constrained  to  a  very 
thin  region  near  the  surface  of  the  airfoil.  For  exam¬ 
ple,  at  Moo  =  0.3,  the  highest  Mach  number  observed 
in  the  contour  lines  of  the  present  experiment  only 
reach  Mi„cai  =  0.84(compared  to  the  Miocai  =  1.3 
inferred  from  surface  pressure  measurements^). 

It  is  also  of  interest  that  this  region  shows  very 
little  longitudinal  gradient,  and  in  fact  shows  the 
constant-density  contours  to  be  almost  parallel  to  the 
surface  of  the  airfoil  in  this  region.  Further  down  the 
airfoil,  these  contours  abruptly  turn  toward  the  sur¬ 
face.  The  shape  of  the  fringe  contours  in  this  region 
suggest  the  presence  of  a  laminar  separation  bubble. 
The  conditions  of  the  present  experiment(NACA  0012 
airfoil,  Re  =  600,000)  are  appropriate  for  a  laminar 
bubble  to  appear  near  the  leading  edge  of  the  airfoil. 
The  formation  of  the  laminar  separation  bubble  was 
also  confirmed  by  the  laser  doppler  velocimetry(LDV) 
studies  performed  in  the  CDSF*^'^^  for  Moo  =  0-3.  k 
=  0  05.  The  angles  of  attack  at  which  the  bubble 
forms  0  =  6®,  and  opens  o  =  14®  are  nearly  the  same 
using  both  LDV  and  PDI.  The  contour  lines  outside 
this  region  curve  more  gradually  toward  the  surface, 
and  can  be  seen  to  define  the  boundary  layer  over  the 
remaining  portion  of  the  airfoil  surface. 

Figure  5  presents  the  flow  field  for  a  =  13.6®.  In 
this  figure,  the  dynamic  stall  vortex  has  formed,  and 
the  density  field  can  be  seen  to  be  dramatically  al¬ 
tered.  The  development  of  the  flow  near  the  surface 
in  this  case  is  of  considerable  interest.  The  most  sig¬ 
nificant  change  in  the  character  of  the  flow  appears 
on  the  upper  surface  near  the  leading  edge.  There 
is  now  a  clearly  defined  shear  layer  leaving  the  lead¬ 
ing  edge;  this  shear  layer  shows  very  strong  gradi¬ 
ents,  as  evidenced  by  the  bending  of  the  light  suffi¬ 
ciently  to  produce  an  opaque  image.  Also,  the  flow 
field  normally  associated  with  the  dynamic  stall  vor¬ 
tex  shows  a  strong  longitudinal  gradient  as  depicted 
by  the  vertical  density  contours  appearing  along  the 
upper  surface  near  25%  chord.  This  figure  presents  a 
quantitative  image  of  the  Mach  number  environment 
of  the  airfoil  at  this  point  in  the  oscillation,  and  can 
be  used  directly  to  evaluate  the  accuracy  of  compu¬ 
tational  models  of  the  dynamic  stall  process. 

4.2  Analysis  of  Stall  Onset  Process 

Detailed  analysis  of  the  interferograms  obtained 
in  the  present  experiment  show  that  the  flow  only 
attains  a  local  Mach  number  of  .83  at  k=.05.  Moo 
=  0.30.  Q  =  10®  -  10®s/na;t.  even  though  in  earlier 
experiments^'  surface  pressures  suggested  local  Mach 
numbers  as  high  as  1.3.  The  two  experiments  are  sim¬ 
ilar  in  many  respects:  both  were  run  at  Moo  =  0.30; 
both  performed  for  q  =  10®  —  10®sinu;t;  both  tests 
were  on  the  NACA0012  airfoil.  However,  the  present 


tests  were  performed  at  a  chord  Reynolds  number 
of  600,000,  while  the  previous  tests  were  performed 
at  Re  =  4x10®.  A  review  of  the  data  from  the  ear¬ 
lier  experiments'*  shows  that  the  NACA0012  airfoil  at 
M=.30,  k=:0.05,  a  =  10®  —  10®smL(;t  reached  a  local 
Mach  number  of  1.22  at  an  a  =  13.9°,  at  which  point 
the  leading  edge  suction  peak  was  lost,  and  the  dy¬ 
namic  stall  process  began.  The  present  experiments 
for  the  same  test  conditions  show  dynamic  stall  to 
begin  at  a  =  12.0®. 

A  laminar  separation  bubble  was  observed  in  the 
earlier  tests'*,  where  it  appeared  as  an  abrupt  change 
in  skin  friction  magnitude  and  frequency  content. 
However,  it  was  only  of  secondary  importance  in  its 
influence  on  the  pressure  distributions  and  appeared 
as  a  small,  localized  “wiggle”  in  the  pressure  traces. 
The  effect  of  the  bubble  was  determined  to  be  sec¬ 
ondary,  since  the  leading  edge  suction  remained  high 
well  after  the  laminar  separation  bubble  was  detected, 
and  the  loss  of  suction  peak  appeared  at  the  leading 
edge,  well  ahead  of  the  laminar  bubble,.  As  noted 
above,  a  laminar  separation  bubble  has  also  been 
identified  in  the  present  set  of  experiments.  However, 
in  the  present  case,  the  bubble  remains  a  significant 
factor  in  the  dynamic  stall  development  -  in  fact,  the 
bubble  behavior  seems  to  be  the  first  indicator  that 
the  dynamic  stall  process  has  begun.  The  character 
of  the  laminar  separation  bubble  changes  significantly 
as  the  dynamic  stall  process  develops.  The  bubble  ap¬ 
pears  at  at  angles  above  o  =  6®  across  a  wide  range  of 
frequencies  and  Mach  numbers.  Although  its  appear¬ 
ance  is  delayed  by  oscillation  at  q  =  8®(see  follow¬ 
ing  section),  it  is  essentially  unaffected  by  frequency 
or  Mach  number  at  any  higher  angle  of  attack.  The 
constant-density  fringes  emanating  from  the  leading 
edge  curve  flatten  significantly  in  this  region,  indicat¬ 
ing  an  alleviation  of  the  local  pressure  gradient.  It 
is  possible  that  since  this  bubble  is  larger  than  in  the 
previous  tesf*,  its  influence  on  the  displacement  thick¬ 
ness  of  the  boundary  layer  is  the  source  of  the  drop 
in  maximum  Mach  number  that  has  been  observed 
throughout  this  set  of  experiments. 

The  bubble  remains  intact  up  to  the  dynamic  stall 
angle.  Figure  6  presents  an  interferogram  illustrat¬ 
ing  the  typical  bubble  geometry,  for  Moo  =  0.33,  k 
=  0.075,  and  a  =  9®(below  the  static  stall  angle). 
This  figure  shows  the  fringes  near  the  airfoil  are  par¬ 
allel  to  the  surface,  followed  by  a  sharp  curving  to¬ 
ward  the  surface  at  the  end  of  the  bubble.  The  pres¬ 
sure  distribution  determined  from  this  interferogram 
is  shown  in  Figure  7.  The  maximum  suction  seen 
in  this  figure  is  approximately  -2.2,  and  a  plateau 
appears  following  the  suction  peak,  extending  from 
.015  <  x/c  <  0.05(the  length  of  the  bubble).  The 
maximum  Cp  observed  in  this  figure  is  significantly 
less  that  that  suggested  by  surface-mounted  instan¬ 
taneous  pressure  transducers  due  to  lack  of  detail  in 
this  region  of  the  photograph;  further  research  will 
focus  on  this  region  in  order  to  quantify  the  details  of 
the  suction  peak.  A  region  of  adverse  pressure  gra¬ 
dient  develops  for  0.05  <  x/c  <  0.08,  in  which  the 
dCp/dx  values  reach  a  high  of  250  while  the  Cp  in¬ 
creases  to  approx  -1.5.  The  pressure  then  gradually 
increases  over  the  rest  of  the  airfoil,  in  a  pattern  sim¬ 
ilar  to  that  of  steady  flow.  For  comparison,  an  inter¬ 
ferogram  obtained  at  a  =  12.4®  is  presented  in  Figure 
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8.  Here,  the  fringe  pattern  is  similar  near  the  lead¬ 
ing  edge,  and  the  beginning  of  the  bubble,  but  the 
downstream  edge  of  the  bubble  now  has  fringes  ex¬ 
tending  downstream,  parallel  to  the  surface,  in  con¬ 
trast  to  the  sharp  turning  toward  the  surface  seen 
in  Figure  6.  In  fact,  the  first  indication  of  break¬ 
down  of  the  flow  leading  to  dynamic  stall  appears  as 
this  change  of  the  fringe  contours  at  the  downstream 
end  of  the  bubble,  suggesting  a  failure  of  the  tur¬ 
bulent  flow  to  reattach.  Finally,  Figure  9,  at  a  = 
13. 07®,  shows  the  beginning  of  the  dynamic  stall  vor¬ 
tex,  which  develops  from  the  extended  bubble  shown 
in  Figure  8.  The  fringes  near  the  surface  downstream 
of  the  bubble  have  now  turned  sharply  toward  the 
surface.  Since  each  fringe  indicates  a  specific  density, 
this  image  shows  that  there  is  now  a  region  of  low 
density,  and  therefore  pressure,  extending  along  the 
airfoils  surface.  This  also  affects  the  flow  away  from 
the  surface,  as  can  be  seen  in  the  outer  fringes  which 
also  curve  around  this  low  pressure  region.  This  pat¬ 
tern  of  opening  of  the  laminar  separation  bubble  was 
detected  throughout  the  frequency  and  Mach  num¬ 
ber  range  explored  to  date,  and  suggests  that  the  dy¬ 
namic  stall  process  in  the  present  experiment  is  either 
caused  by  a  breakdown  of  the  laminar  separation  bub¬ 
ble,  or  by  events  occurring  within  the  bubble  which 
have  yet  to  be  documented.  This  is  of  significant  in¬ 
terest,  since  it  re-emphasizes  that  dynamic  stall  can 
be  caused  by  a  wide  range  of  fluid  physics  while  result¬ 
ing  in  very  similar  global  flow  behavior.  Future  tests 
in  the  present  program  will  investigate  various  types 
of  leading-edge  boundary  layer  transition  generators, 
so  that  the  CDSF  airfoil  can  fully  reproduce  the  var¬ 
ious  types  of  leading-  and  trailing-edge  dynamic  stall 
conditions  observed  in  full  scale  tests. 

4.3  Influence  of  Reduced  Frequency  and  Mach  Number 

REDUCED  FREQUENCY:The  effect  of  reduced 
frequency  on  the  flow  developing  on  the  airfoil  can  be 
seen  in  Figure  10.  In  this  figure,  the  influence  of  un¬ 
steadiness  is  shown  at  each  of  four  angles  of  attack. 
The  top  row  of  photographs  shows  the  effect  of  fre¬ 
quency  on  the  flow  at  q  =  8®.  The  presence  of  the 
laminar  separation  bubble  can  be  seen  at  k  =  0.0. 
At  this  angle  of  attack,  the  bubble  is  suppressed  by 
oscillation;  the  flow  at  k  =  0.075  and  0.10  do  not 
show  the  bubble  -  the  contours  are  similar  to  those  at 
lower  angle  of  attack(not  shown).  At  a  =  10®(second 
row  of  photographs)  the  bubble  is  present  throughout 
the  frequency  range.  At  a  =  12®(third  row  of  pho¬ 
tographs),  the  first  indication  that  oscillation  delays 
stall  is  observed:  the  steady  flow  is  stalled,  but  even 
low  level  oscillation(k  =  0.025)  suppressed  separation 
and  re-establishes  fully  attached  flow.  At  q  =  14®, 
the  influence  is  even  more  clearly  discernable:  as  the 
frequency  increases  to  k  =  0.05  and  higher,  the  stall  is 
suppressed,  and  a  dynamic  stall  vortex  appears.  The 
development  of  this  vortex  is  delayed  as  the  frequency 
is  increased.  Thus,  as  the  frequency  reaches  k  =  0.10, 
the  dynamic  stall  process  is  just  beginning  at  an  angle 
two  degrees  higher  than  the  angle  at  which  the  airfoil 
was  completely  stalled  in  steady  flow. 

MACH  NUMBER:  The  effect  of  Mach  number  at 
k  =  0  05  can  be  seen  in  Figure  11.  At  q  =  6®(first 


row  of  photographs),  the  flow  varies  smoothly  from 
the  stagnation  point  to  the  trailing  edge(out  of  the 
view  of  the  photograph).  As  the  Mach  number  is  in¬ 
creased,  a  progressive  increase  in  the  fringe  count  is 
apparent,  as  can  be  expected,  due  to  the  greater  den¬ 
sity  variations  that  occur  at  the  higher  Mach  num¬ 
bers.  At  a  =  8®  (second  set  of  photographs),  the 
presence  of  the  laminar  separation  bubble  is  observed. 
This  bubble  is  not  measurably  altered  as  the  Mach 
number  is  increased,  except  for  the  increase  in  fringe 
density  reflecting  the  increase  in  Mach  number.  At  a 
=  10®(third  row^  some  very  interesting  flow  physics 
can  be  seen  in  the  laminar  separation  bubble.  As  the 
Mach  number  increases  to  Moo  =  0.42,  vertical  lines 
appear  in  the  bubble.  These  vertical  lines  appear  even 
more  clearly  at  Moo  =  0.47,  where  the  fringes  are  seen 
to  bend  discontinuously  as  they  pass  through  the  ver¬ 
tical  lines.  The  local  Mach  number  obtained  by  fringe 
counting  was  found  to  be  1.0  -  1.1;  thus,  the  local  flow 
is  slightly  supersonic.  The  multiple  shock  waves  seen 
are  a  feature  of  the  transonic  flow  over  mildly  curving 
surfaces;  this  condition  of  multiple  shocks  have  been 
observed  in  other  experiments^®' It  is  very  interest¬ 
ing  that  the  current  study  shows  the  phenomenon  to 
be  present  even  in  unsteady  flows,  where  the  events 
occur  at  different  angles  of  attack  depending  on  the 
reduced  frequency. 

At  a  =  12®(fourth  row),  the  influence  of  com¬ 
pressibility  on  the  dynamic  stall  process  is  quite  evi¬ 
dent.  The  flow  on  the  airfoil  at  Moo  —  0.28  is  clearly 
attached,  as  is  the  flow  at  Moo  =  0.30.  However,  at 
Moo  =  0.33,  the  dynamic  stall  process  has  begun,  and 
at  Moo  =  0.38  the  dynamic  stall  vortex  is  evident.  At 
Moo  =  0.42,  the  vortex  is  already  past  x/c  =  .25,  and 
at  Moo  =  0.47,  the  vortex  is  already  out  of  the  field  of 
view.  These  photographs  graphically  show  that  the 
dynamic  stall  inception  angle  is  a  direct  function  of 
the  Mach  number. 

5.  QUANTITATIVE  ANALYSIS  OF 
THE  UNSTEADY  FLOW  FIELD 

5.1.  Mach  Contours 

The  interferograms  presented  in  this  study  offer 
quantitative  information  concerning  the  Mach  con¬ 
tours  of  the  flow  under  investigation. 

For  a  standard  interferometer,  the  path  length 
difference  APL  due  to  density  (or  phase)  changes  can 
be  related  to  the  fringe  number  e  as, 

where  n  is  the  refractive  index  in  the  signal  beam, 
Uref  is  that  in  the  reference  beam  and  Aq  is  the  wave 
length  of  the  light  used. 

For  a  two-dimensional  flow,  the  above  equation 
can  be  simplified  to, 

f  =  (n-  nr«/)-r- 

where  L  is  the  test  section  span.  If  f  is  an  integer, 
then  the  fringe  is  bright  and  if  it  is  a  half  integer,  the 
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fringe  is  dark.  Using  the  Gladstone  -  Dale  equation 
and  the  perfect  g^ls  equation,  this  equation  can  be 
reduced  to  (See  Goldstein,  Reference  15) 

P-PreS=i-^)C{)  =  Ae 

no  —  1  L 

Since  Ac  =  532  nm,  L  =  25cm,  (no  —  1)  =  2.733 
X  10“'*  (from  Goldstein),  and  po  =  l.2lkg/Tn^,  the 
constant  A  can  be  determined.  For  the  specific  case 
of  the  present  experiments, 

p  —  Po  =  0.009421e 
or 

£.  =  +  ^ 

Po  Po  po 

Since  ^  is  a  function  of  the  free  stream  Mach 

fio 

number  only,  ^  can  be  determined  by  knowing  the 
fringe  number.  Note  that  in  this  method,  a  positive 
fringe  number  represents  deceleration  and  vice  versa. 

5.2  Digitization  of  Interferograms 

In  order  to  make  full  use  of  the  quantification 
that  interferometry  offers,  the  photographs  contain¬ 
ing  the  interferograms  were  digitally  processed,  and 
the  leading  edge  region  was  greatly  enlarged.  Three 
registration  markers  were  placed  on  one  of  the  rotat¬ 
ing  tunnel  windows  which  support  the  airfoil  so  that  a 
proper  orientation  and  scaling  of  the  fringe  patterns 
could  be  accomplished( Figure  12).  The  two  trian¬ 
gles  below  the  airfoils  are  arranged  so  that  the  dis¬ 
tance  between  the  vertical  faces  is  exactly  20%  chord; 
the  horizontal  surfaces  of  these  triangles  are  placed 
so  that  a  line  connecting  these  surfaces  is  parallel  to 
the  airfoil  chord,  and  10%  below  the  center  line  of  the 
airfoil.  The  triangle  above  the  airfoil  is  placed  so  that 
its  lower  surface  is  20%  chord  from  the  top  surface  of 
the  triangle  located  below  the  airfoil.  A  line  connect¬ 
ing  the  vertical  surfaces  of  these  two  triangles  will  be 
perpendicular  to  the  chord  at  the  25%  chord  location. 

Each  interferogram  was  digitally  scanned,  and 
scaled  to  a  selected  size  by  use  of  these  registration 
•markers.  A  grid  was  superimposed  on  the  digitized 
image,  and  the  grid  block  containing  the  leading  edge 
region  was  selected  for  further  enlargement.  This  im¬ 
age  was  then  processed  using  an  adaptive  contrast- 
enhancement  procedure  specially  developed  for  this 
purpose.  This  procedure  subdivided  the  region  be¬ 
fore  applying  the  enhancement,  and  then  locally  opti¬ 
mized  the  contrast.  The  locally  optimized  contrast  in 
each  subregion  was  then  smoothly  matched  to  the  ad¬ 
joining  regions,  resulting  in  a  greatly  improved  inter¬ 
ferogram  for  digital  fringe  mapping( Figure  13).  This 
image  shows  a  clear  delineation  of  the  airfoil  surface. 
This  clarity  was  obtained  by  using  a  photograph  of 
the  airfoil  taken  under  no-flow  conditions,  with  the 
PDI  plate  removed.  This  no- flow  airfoil  image  was 
then  digitized,  scaled,  and  superimposed  on  the  dig¬ 
itized  image  of  the  interferogram,  and  the  combined 
digital  image  was  then  used  for  further  analysis. 

The  combined  image  was  then  mapped  using  a 
specially-designed,  screen-oriented  digitizing  program 


which  interactively  permitted  tracing  of  the  various 
fringes  by  movement  of  the  cursor  via  a  computer 
“mouse”.  Samples  of  the  resultant  fringe  contour 
maps  can  be  seen  in  Figures  14  and  15. 

5.3  Instantaneous  Pressure  Distributions 

Each  fringe  in  an  interferogram  is  a  line  of  con¬ 
stant  density,  and  by  use  the  procedure  outlined  in 
Section  5.1,  denotes  a  constant  Mach  number.  There¬ 
fore,  the  fringe  maps  described  above  can  be  used  to 
quantify  the  effect  of  unsteadiness  on  the  local  pres¬ 
sure  distributions  of  the  oscillating  airfoil  studied  in 
the  present  experiments.  Figure  16  presents  pressure 
distributions  near  the  leading  edge  of  the  airfoil  for 
four  steady  angles  of  attack  at  A/oo  =  0.33.  These 
plots  are  obtained  by  mapping  the  x-location  of  the 
intersection  of  each  fringe  with  the  airfoil  surface.  The 
fringe  numbers  are  then  converted  to  Cp,  and  the 
pressure  distributions  calculated.  In  each  case  shown, 
the  distributions  are  mapped  up  to  the  point  where 
the  last  fringe  intersects  the  surface  before  the  laminar 
separation  bubble  occurs,  and  thus  show  the  pressure 
distribution  associated  with  the  laminar  bubble  for- 
mation(except  for  a  =  4°,  where  the  bubble  was  not 
present). 

Figure  17  presents  the  pressure  distributions  at 
k=0.05  and  k=0.10  at  a  =  S'*.  The  pressure  distri¬ 
butions  for  steady  flow  at  o  =  6®  and  8®  are  included 
for  comparison.  TThe  curves  for  k  =  0.05,  and  0.10  are 
flatter  and  thus  show  progressively  lower  gradients. 
The  pressure  distribution  for  k  =  0.10,  at  o  =  8®  is 
virtually  identical  to  the  distribution  for  k  =  0.0,  o 
=  6®,  showing  that  unsteadiness  has  indeed  reduced 
the  instantaneous  local  pressure  gradients.  This  ex¬ 
perimental  evidence  of  the  lessening  of  local  adverse 
pressure  gradient  is  the  first  data  to  show  this  detail  of 
the  developing  pressure  gradients,  since  the  cost  and 
complexity  of  installation  of  instantaneous  pressure 
transducers  in  sufficient  number  to  show  this  detail 
has  precluded  such  endeavors  in  the  past. 

Figure  18  presents  the  maximum(suction)  Cp  val¬ 
ues  as  a  function  of  the  angle  of  attack  at  Moo  = 
0.33,  k=  0.075.  The  maximum  Cp  was  determined 
by  counting  the  fringes  through  the  stall  vortex  and 
tracing  them  back  to  the  leading  edge;  in  the  attached 
flow  near  the  leading  edge,  any  particular  fringe  can 
be  assigned  a  Cp  value  using  isentropic  flow  relations. 
This  plot  shows  that  suction  steadily  develops  over 
the  airfoil  as  the  angle  of  attack  increases.  At  a  = 
12.5®,  the  Cp  has  reached  a  value  =  -4.7,  and  remains 
at  this  level  till  q  =  13.2®,  where  it  starts  decreasing. 
It  is  believed  that  the  dynamic  stall  process  begins  at 
a  =  12.5°  in  this  case;  during  the  vortex  formation 
the  suction  peak  remains  constant.  However,  once  the 
vortex  begins  to  convect,  the  suction  peak  drops.  It  is 
interesting  that  even  though  the  suction  peak  drops, 
there  is  a  substantial  suction  present  over  the  airfoil 
as  the  vortex  grows,  and  the  increased  extent  of  the 
low  pressure  region  enables  generation  of  higher  lift 
until  the  vorticity  leaves  the  airfoil. 

LIMIT  ON  Cp  MIN:  Also  plotted  in  Figure  18 
are  the  minimum  Cp  values  for  Moo  =  0.3  at  k=0.05, 
0.075,  and  0.10,  when  they  are  first  observed(i.e.  be¬ 
fore  the  vortex  could  be  detected  in  the  interfero- 
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grams).  It  is  seen  that  Cpmin  =  -4.96  for  all  three 
cases  at  Moo—  0.3.  The  most  noteworthy  difference 
is  that  2is  the  reduced  frequency  increases,  the  an¬ 
gle  of  attack  at  which  the  peak  suction  is  encoun¬ 
tered  and  the  stall  process  begins  also  incre2ises.  For 
example,  a,taii  =  12.13°  ,  12.5°,  and  13.14°  for  k 
=  0.05,  0.075,  and  0.10  respectively.  These  results 
clearly  suggest  that,  at  a  given  Mach  number,  stall 
develops  once  the  Cp  value  reaches  a  unique  value, 
independent  of  reduced  frequency.  Since  the  value  is 
reached  at  progressively  higher  angles  with  increasing 
frequency,  oscillating  an  airfoil  at  higher  frequencies 
delays  the  stall  to  higher  angles.  As  has  been  shown 
in  Figure  17,  the  pressure  gradients  are  also  reduced 
by  oscillation,  resulting  in  lowering  of  the  equivalent 
instantaneous  angle  of  attack.  Thus,  at  least  part  of 
the  delay  of  stall  can  be  viewed  as  a  modification  of 
the  inviscid  flow  imposed  on  the  boundary  layer.  The 
effect  of  compressibility  on  this  flow  development  is 
apparent  in  that  as  the  Mach  number  is  increased, 
the  Cpmin  value  decreases,  as  can  be  seen  in  the  fig¬ 
ure  for  Moo  —  0.4,  k  =  0.05,  where  a  value  of  -3.9  is 
attained  at  q  =  11.2°. 

6.  CONCLUDING  REMARKS 

Point  diffraction  interferometry  has  permitted  de¬ 
tailed  study  of  the  complex  unsteady  flow  near  the 
leading  edge  of  an  oscillating  airfoil,  and  quantitative 
flow  information  has  been  obtained  both  on  the  sur¬ 
face  and  in  the  surrounding  flow  field  for  a  range  of 
frequencies  and  Mach  numbers.  A  laminar  separation 
bubble  is  observed  in  most  of  the  higher  angle  con¬ 
ditions,  although  the  occurrence  of  the  bubble  can 
be  delayed  by  unsteadiness.  Locally  supersonic  flow 
was  observed  near  the  leading  edge,  but  the  region 
of  supersonic  flow  is  quite  small.  Digitization  of  the 
interferograms  has  permitted  detailed  mapping  of  the 
instantaneous  pressure  distribution  near  the  leading 
edge,  and  has  shown  that  unsteadiness  significantly 
relieves  the  pressure  gradients  that  occur  in  this  re¬ 
gion. 
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Fig.  3.  Sdipmatic  of  thr  Layout  of  Optics  for  Point-Diffraction  Interferometry. 
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Figure  1 1(  eonrliuled  )  Point  Dii[r;ii  ti<iii  Iiiterferogrniiis  Sliowing  Effect  of  Mach  N\imt>er 
on  Dynamic  Stall  Develojuuent  at  Four  .Angles  of  .Attack.  ;it  k  --  O.O-j. 
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Figure  12.  Diagram  of  Registration  Markers  placed 
on  CDSF  Tunnel  Window. 


Figure  13.  Enlarged  Image  of  Leading  kMge  SelMted 
FVom  InterferogramfAfter  Digital  Enhancement). 


Figure  14.  Fringe  Contours  Near  the  Leading  Edge  of 
an  NACA  0012  Airfoil  at  Af*  =  0.33,  a  =  8.0°,  (a) 
k  =  0.0;  (b)  k  =  0.10. 
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Figure  15.  FVioge  Contours  Neat  the  Leading  Edge  of 
an  NACA  0012  Airfoil  at  Af»  =  0.33,  o  =  10.0°,  (a) 
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Figure  16-  Pressure  Distributions  for  Several  Angles 
of  Attack  for  NACA  0012  Airfoil  in  Steady  Flow  as 
Determined  From  PDl  Interferograms,  Afoo  =  0.30,  k 
=  0.0. 


Figure  17.  Pressure  Distributions  for  Several  Reduced 
Frequencies  for  NACA  0012  Airfoil  as  Determined 
From  PDI  Interferograms,  Moo  =  0.30,  Q= 
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Figure  18.  Minimum  Cp  as  Function  of  Angle  of  At¬ 
tack  at  Various  Mach  Numbers,  as  Determined  From 
PDI  Interferograms. 


